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ABSTRACT 

We present an empirical determination of the white dwarf cooling sequence in the globular cluster 47 
Tucanae. Using spectral models, we determine temperatures for 887 objects from Wide Field Camera 
3 data, as well as 292 objects from data taken with the Advanced Camera for Surveys. We make the 
assumption that the rate of white dwarf formation in the cluster is constant. Stellar evolution models 
are then used to determine the rate at which objects are leaving the main sequence, which must be 
the same as the rate at which objects are arriving on the white dwarf sequence in our field. The result 
is an empirically derived relation between temperature (T e //) and time (t) on the white dwarf cooling 
sequence. Comparing this result to theoretical cooling models, we find general agreement with the 
expected slopes between 20,000K and 30,000K and between 6,000K and 20,000K, but the transition 
to the Mestel cooling rate of T e ff oc t~ 0A is found to occur at hotter temperatures, and more abruptly 
than is predicted by any of these models. 



1. INTRODUCTION 

The majority of stars that form in our Galaxy will 
eventually end up as white dwarfs. Given enough mass 
to fuse hydrogen but not enough mass to eventually 
form a neutron star or a black hole, the range of ini- 
tial masses for white dwarf stars encompasses much of 
the lower end of the initial mass function. The study 
of such objects can therefore reveal a wealth of informa- 
tion about the population that formed them. A predic- 
tive model of white dwarf cooling can als o be used to 
deriv e precise ages for ste llar populations (| Oswalt et al.l 
119961: lHansen et all l2007h . Additionally, white dwarfs 
themselves can serve as laboratories for constraining the 
physics that determine their cooling rates. 

The aim of this work is to obtain an empirical mea- 
surement of the rate at which white dwarfs cool. Theo- 
retical work in this field began with Mcstcl's 1952 paper 
in which he layed out a simple model for the rate at 
which white dwarfs should lose energy. This model as- 
sumes a completely degenerate and isothermal core sur- 
rounded by a thin non-degenerate envelope. Thermal 
motion of the non-degenerate ions as well as the slow 
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collapse of the star itself releases gravothermal energy 
from the core of the white dwarf, while the rate of en- 
ergy loss fro m the star is determined by the opacity in 
the envelope (|Mestel|[l952l) . The end result is a predic- 
tion that the temperature of a white dwarf should de- 
crease with time since formation, t, as t~ 0A . This basic 
model has been gradually improved upon since then, with 
the introduction of more and more physical considera- 
tions. Summaries of progress i n this field can be found in 
iD'Antona fc MazzitelU] (1199011 Hansen! (|2004ft . and most 
recently lAlthaus et al.l ( 201C ) . Hansen also presents a 
comparison of some modern cooling models. However, 
comparisons of white dwarf cooling models to observable 
populations tend to focus on the cool end of the sequence, 
which is the most relevant for dating. We present here 
a detailed comparison of a number of current models to 
measurable cooling rates up to temperatures of 30,000K. 

In this paper, we will analyze a group of white dwarfs 
that have all evolved from the same initial population. 
These objects are part of the Galactic globular cluster 
47 Tucanae (NGC 104). All of the objects in this cluster 
formed at the same time, or at least formed over a period 
of time that is very small when compared with the cur- 
rent age of the cluster. The fact that our entire sample 
of white dwarfs comes from a single population allows us 
to draw a number of simplifying assumptions. 

The first is that all of the white dwarfs are of a similar 
mass. Because all of the stars should have formed around 
the same time, all of the objects on the main sequence 
have had roughly the same amount of time to evolve. 
This implies that the variance in the masses of objects 
currently leaving the ma in sequence should be small. Ini- 
tial final mass relations (jKalirai et alll2009| ) then suggest 
that variance in the mass of the white dwarfs, which these 
turn-off stars will eventually become, is also small. This 
agrees with the analysis of white dwarf spectra in glob- 
ular clusters, which suggest that all of the objects are o f 
similar mass (jKalirai et all 120091 : iMoehler et all l2004f l . 
It is important to note that this assumption will not be 
true for very old white dwarfs, but since our entire sam- 
ple should be younger than a few Gyrs, this assumption 
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is valid. 

The second assumption is that objects arrive on the 
white dwarf sequence at a constant rate. This follows 
from the fact that the rate at which objects are leaving 
the main sequence is constant. This is supported by a 
comparison of the luminosity function of the lower giant 
branch with a model luminosity function with a constant 
evolution rate (see Section [5J . Given that this rate is 
constant, and all objects leaving the main sequence will 
eventually become white dwarfs, the rate at which ob- 
jects are arriving onto the white dwarf sequence must 
also be constant. 

These two points allow us to put all of our white dwarfs 
on the same cooling sequence by fitting models to each 
object's spectral energy distribution (SED) to determine 
a temperature, and then requiring that the white dwarfs 
are uniformly distributed in time along the sequence. 
This approach allows us to measure the rate at which 
the temperature of the white dwarfs is changing in a way 
that is independent of any white dwarf cooling models. 

2. DATA 

Our data set consists of 121 orbits of Hubble Space 
Telescope (HST) observations. During these orbits, ex- 
posures were taken simultaneously using the Advanced 
Camera for Surveys (ACS) and Wide Field Camera 3 
(WFC3). The current paper focuses on the WFC3 UVIS 
data in F390W and F606W. These data are comprised 
of 26 orbits covering 13 adjacent fields (1 orbit in each 
filter). The description of the methodology in this paper 
will refer to that data set. Similar methods were used 
to produce a cooling curve from the ACS data, which is 
taken in F606W and F8UW. The WFC3 data cover 
over 60 arcmin 2 over a range of 6.5-17.9 pc from the core 
of 47 Tuc. The ob servations and data reduction are de- 
scribed in detail in IKalirai et al.l (|2012D . 

3. SELECTING OUR WHITE DWARF SAMPLE 

Selection began with a catalog g enerated from the fi - 
nal drizzled images as described in IKalirai et al.l (|2012f ) . 
RMS residua ls of the PSF fit from DAOPHOT I I 
(lStetsonlll987ri . SHARP from ALLSTAR (lStetsonlll994r) . 
and photometric errors from DAOPHOT were used to 
cut down the full catalog to a higher quality sample. 
Stars that made the cut were required to have errors 
of less than 0.1 magnitudes in ^3901^ and F606W, a 
SHARP value within 2a of the mean SHARP value, and 
an RMS of the PSF fit less than 2.5%. The CMD after 
selection cuts are made is shown in Figure [T] 

The region around the white dwarf sequence be- 
tween roughly 22nd and 27th magnitudes in F606W was 
bounded manually. The objects inside this region were 
used as our final sample. Although some real objects 
are almost certainly lost due to these cuts, this can be 
properly accounted for during the analysis of artificial 
star tests. The same cuts are used to determine which 
objects are "found" from our input artificial stars. 

Uncertainty also arises from manually bounding our 
sequence. We analyze the distribution in color along the 
white dwarf sequence and model the spread in color as 
Gaussian. If we integrate this Gaussian over our enclosed 
area and compare this to the total integral, we find that 
we expect to lose 14 objects due to this bounding region. 
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Fig. 1. — The CMD after making selection cuts in error, sharp- 
ness, and x 2 of the PSF fit. The final WD selection region is 
shown in the enclosed area (red). 887 objects are found within this 
enclosed region. 



This means that of our final sample of 887, each object 
should count as (14+887)/887 objects (prior to incom- 
pleteness corrections). The uncertainty introduced by 
this is of order 1.5% in the time scaling of our object list. 
This scaling is discussed in more detail in Section [5j 

4. FITTING WHITE DWARF TEMPERATURES 

To determine an effective surface temperature (T e ff) 
for e ach object in ou r samp le, we utilize spectral models 
frompTrcmblav et al.l (|2011f) . The models are for thick at- 
mosphere (the H envelope mass fraction is approximately 
c[h = 10~ 4 ) DA white dwarfs, and are parametrized over 
a grid of surface gravity and T e ff values ranging from 
6.0 to 10.0 and 6,000 K to 120,000 K, respectively. The 
DB fraction is assumed to be negligible , whic h is sup- 
ported by analysis from iWoodlev et al.l (|2012f) . These 
models (F moc i) are then combined with the filter through- 
put cu rves (S\), which are available online from[STScI 
(2009), to generate a magnitude in a particular filter for a 
given set of model parameters. This is shown in equation 
ffl 



M m od = -2.5 log! 
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The model magnitudes are fit to our sample of 887 
white dwarfs using a maximum likelihood fitting method. 
The likelihood of the data given the model is: 



Lidata \ parameters) = J^J exp 
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Here the model parameters are (m — M)o, E(B — V), 
Mwd, and T e ff. To condense this to the likelihood of 
T e ff given our data, we integrate L(data \ parameters) 
over (m— M)o, E(B—V), and Mwd with Gaussian priors 
of 13.3 ± 0.15, 0.04 ± 0.02 and 0.53 ± 0.02 M respec- 
tively (jWoodlev et alJ 12012). The dista nce prior used 
here i s the mean of all values quoted in iWoodlev et al.l 
()2012|) that do not depend on the white dwarf sequence. 
The maximum likelihood in T e f f is found for each object. 
The uncertainty in the distance is tw ice as large as the 
value quoted bv IWoodlev et all (|2012l ). This more accu- 
rately reflects the distribution of distance measurements 
made by various groups, as presented in that paper. The 
(Tj in Equation [2] is the estimated photometric error of 
each data point. These values are de r ived f rom artificial 
star tests presented in iKalirai et al.l (|2012l ) . The error 
distributions of the two WFC3 filters are shown in Fig- 
ure El 

5. DETERMINING THE TIMESCALE 

5.1. A Sorted List of Temperatures 

After an effective surface temperature is determined 
from fitting spectral models to the photometry of each 
object, this list must be sorted from hottest to coolest. 
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Fig. 2. — A model WD spectrum with an effective surface temper- 
ature of 10,000K and a surface gravity of log(g)=7.5. Throughput 
curves for our two filters are shown overplotted on a separate axis 
scale. 
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Fig. 3. — The standard deviation of input - output magnitude 
derived from artificial star tests described in Kalirai et al. (2012) 
(for all fields combined). 



This sorted list is the first step toward determining the 
rate at which the white dwarfs are cooling. Because the 
objects are arriving on the white dwarf sequence at a 
constant rate, the position in the sorted list should be 
proportional to the time spent on the sequence. That 
is to say, the time between when the 100th and 200th 
objects arrive on the sequence should be the same as the 
time between the 200th and 300th. 

Completeness of the sample also needs to be taken into 
account. The cooler objects are also fainter and thus less 
likely to be seen in a crowded field. Usin g artificial star 
tests as described in IKalirai et all (|2012l ) , we can assign 
a completeness correction to each object by first binning 
the artificial star data by magnitude and position. The 
completeness correction or fraction in some bin is then: 
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This can be interpolated to the magnitude of each ob- 
ject in our sample and the distance from the cluster cen- 
ter. To calculate the position of each object in a com- 
pleteness corrected sorted list we use the following equa- 
tion: 



x(n) 



i=0 



(i). 



(4) 



It is easy to see that if the entire sample is 100% com- 
plete, then this is just an array of positions in the sorted 
list (1, 2, 3, etc). In our sample the completeness correc- 
tions are small. Even the faintest objects are more than 
75% complete, and overall the sample is 89% complete. 
The completeness fraction as a function of magnitude is 
shown in Figure |U The artificial star tests were actu- 
ally analyzed field by field. Fields closer to the center 
of the cluster are obviously more dense, and have larger 
incompleteness. The results in Figure @] are the total for 
all fields. Incompleteness was also considered in both 
filters. It was found to make no difference whether in- 
completenes values from a single filter, or a combination 
of both filters was used. Since the corrections were so 
small, this had no effect on our final fitted parameters. 
Completeness corrections for the ACS data were calcu- 
lated in the same way. The entire ACS sample is greater 
than 90% complete. 

A sorted, completeness corrected list is shown in log- 
log space in Figure [5] This relation is well fit with a 
broken power law. The slopes and "break" of such a 
fit are also shown in the Figure. These parameters do 
not depend at all on the absolute scaling of the x-axis, 
since this only acts as an additive constant in log-space 
(logio [At/ object]). These power law parameters and 
their uncertainties were derived by iterating our tem- 
perature fitting analysis. In each iteration we assume 
parameters of (m-M)o, E(B-V), and Mwd with values 
drawn from the prior distributions discussed above. Af- 
ter each iteration, the temperatures are sorted and fit 
with this broken power law. The distributions of these 
fit parameters define the best fitting value, and the uncer- 
tainties. This process was also repeated for the ACS field 
with 292 WDs. The results are shown in the lower panel 
of Figure [5] All fitted parameters are consistent with the 
parameters from the WFC3 data. The uncertainties in 
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Fig. 4. — Completeness fraction of our WD sample as a function 
of F606W magnitude (averaged over all fields). 
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Fig. 5. — The fitted T e ff as a function of object number for our 
WDs. The object number is proportional to the time on the white 
dwarf sequence, and is calculated as shown in equation [4] The y- 
axis is just the temperature fit to each object's SED, as discussed 
in section [J] The results from both data sets agree well within 
uncertainties. Residuals are shown as Dato j - I f° 
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these fit parameters take into account the uncertainty in 
the photometry of our individual objects, as well as the 
uncertainties in all of our marginalized parameters. 

5.2. The Rate at Which Stars are Leaving the Main 
Sequence 

In order to put the completeness corrected list into 
units of time, so that our cooling curve can be com- 
pared with theoretical models in both dimensions, we 
need to find the rate at which objects are arriving onto 
the white dwarf sequence in our field. This rate should 
be the same as the rate at which objects are leaving the 
main sequence. One could imagine that dynamical re- 
laxation of white dwarfs could cause problems with this 
assumption. The white dwarfs are considerably less mas- 
sive than their giant branch progenitors and should be 
expected to relax on timescales similar to the relaxation 
time for the field. However, this relaxation time is long. 
At the half-li ght radius the relaxation time is 3.5 Gyrs 
(|Harrislll996l ). Usin g the velocity dis persion profile of 
47 Tuc calculated bv lLane et al.l (|2010() . and the density 
profile from our own data, we calculate an expected mean 
relaxation time for our WFC3 field to be 10 Gyrs. We 
can compare this to the expected time since the white 
dwarfs in our sample were the mass of current turn-off 
stars, which is < 2.5 Gyrs for st ars brighter tha n 27.8 in 
^390^ (utilizing models from iFontaine et al.l (|2001ft ). 
Given that these objects have only existed for a small 
fraction of a relaxation time at their current mass, there 
should be no observable mass segregation. Additionally, 
comparing the radial distributions of the white dwarfs 
and the giant branch stars, no significant difference is 
apparent, as shown in Figure [6] This indicates that the 
white dwarfs have not yet relaxed, and their distribution 
in the field should be the same as their progenitors. To 
determine the age scaling for our cooling curve we need 
to analyze the giant branch of the CMD. The upper gi- 
ant branch i s saturated in our exposures, so we use the 
catalog from iStetsonl ([2000) of an overlapping field in B 
and V, selecting only those objects that would be found 
within the region of our WFC3 observations. 
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Fig. 6. — A comparison of the cumulative radial distributions of 
the WD and lower RGB stars. A KS test indicates that there is 
a 73% chance that two samples of this size drawn from the same 
underlying distribution would differ by at least this much. We 
therefore find no evidence of significantly different radial distribu- 
tions for the WD and RGB stars. 

On this portion of the CMD, stars evolve almost verti- 
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cally as their outer layers expand accompanied by a very 
moderate decrease in effective surface temperature. To 
determine an evolution rate from these obj ects, we use 
theoretical models from lDotter et ahl (pOOl . Our input 
parameters for the model are [Fe/H], (m-M)o, current 
cluster age (^4), and an age range {AA). This age range 
corresponds to the time it takes an object to evolve from 
one point on the lower giant branch to another point 
(point Pi to point P2 in Figurc[7]). In our case, all of these 
parameters except the age range are well determined. We 
calculate the likelihood of our magnitude range resulting 
from this model over this whole parameter space using 
Equation [SJ We marginalize out all parameters except 
the age range, using priors that describe the measured 
values of [Fe/H], (m-M) , and A. 
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AM r is the selected magnitude range (shown in Fig- 
ure [7]), while AM m is the range predicted by the 
model given a part i cular choice of input parameters. 
IChabover fc Kraus s (2002) show that there is an un- 
certainty of order 3% in the age of such models. To 
first order, object age and magnitude are linearly re- 
lated in this region, and so this translates to an un- 
certainty of 3% in the output magnitude range. This 
uncertainty is used as a in Equation [5] The priors 
used in Equation [5] are all Gaussian with the following 
means and standard deviations: [Fe/H] = —0.75 ±0.08, 
(to - M) F ^ mw = 13.42 ± 15, and A = 11.5 ± 0.75 
Gvr (fM c Willia m fc BernsteirJl200l IWoodlev et aDIMl 
iDotter et al.ll2009D . 

The result is a likelihood distribution of AA for our 
chosen magnitude range that takes into consideration un- 
certainties in the input model parameters, as well as un- 
certainty in the model itself. This distribution is shown 
in Figure [8] 

Our model prediction for the timespan corresponding 
to this chosen magnitude range is 496±lf Myrs (2a un- 
certainties). The uncertainties on this value are entirely 
due to uncertainties in the model as well as the input 
parameters. To determine the rate at which objects are 
evolving off the main sequence, the fit time interval must 
be divided by the number of objects (N = 287) in this 
region of the CMD (shown in Figure [7]) within our field. 

The uncertainty in the number of objects must be y/N, 
resulting in a final rate of 1.7 ± 0.2 Myrs per object, an 
age scaling uncertainty of 12%. This value has no phys- 
ical meaning for the cluster as a whole. It is dependant 
entirely on our field. However, as this is the same field 
as our WD sample, it is appropriate to assume that this 
rate is the same as the rate at which objects are arriving 
on the WD sequence. This time/object is then multiplied 
by our ordered and completeness corrected list, resulting 
in an empirically measured relation between the temper- 
ature of the white dwarfs, and their age on the sequence. 
This method also assumes that the rate of evolution off 
the main sequence is constant. This assumption is well 
supported by analyzing the cumulative distribution of 
magnitudes between the two points in Figure [7] This 
distribution is shown compared against a model distri- 




FlG. 7.— The CMD from [Stetson (2000), including only ob- 
jects that fall within the region of our WFC3 data. The region 
of the CMD betwe en points Pi and P2 is fit with models from 
IDotter et al. (200H) to determine the time it takes an object to 
evolve between the points. 
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Fig. 8. — Likelihood of AA fit to our selected magnitude range, 
after marginalizing out all other model parameters. AA corre- 
sponds to the time it takes an object to evolve from point Pi to 
point P2 in Figure [7] 

bution assuming a constant evolution rate in Figure [9] 

6. THE COOLING CURVE 

The cooling curve determined from our data is shown 
in Figure [TO] The shape of the cooling curve is well con- 
strained, as can be seen in the residuals of Figure [5] The 
uncertainties quoted there represent the la uncertain- 
ties in the upper and lower slopes taking into considera- 
tion both the random uncertainties in the magnitudes of 
each object, as well as the unce rtainty in the magnitude 
calibration (Kalirai et al. 2 0121) . The contours shown in 
Figure [TU] display the systematic uncertainty in the x- 
direction that arises from our uncertain time-scaling fac- 
tor. This is the dominant source of uncertainty in our 
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Fig. 9. — A cumulative distribution of Stetson's RGB magnitudes 
between Pi and P2 in Figure [7] compared to the expected distri- 
bution assuming a constant evolution rate. A KS test suggests 
that there is a 40% chance that a sample drawn from the model 
distribution would differ by at least as much as this data does. 

time scaled cooling curve. 

6.1. The Models 

We com pare our data to mode ls from four groups. The 
models of Fontaine et al.l ()2001l ) have pure C cores with 
H and He envelo pe mass fract i ons o f qn = 10~ 4 and 
qHe = 1CT 2 . The (Hansen et al.l ()2007l ) models also have 
envelope mass fractions of qn = 10~ 4 and qHe = 10~ 2 , 
and a mixed C/O cor e with profiles from | Salaris et al.l 
(|1997D . Models from lLawlor fc MacDonald! ()2006D are 
evolved from the main sequence. By the time they reach 
the white dwarf stage, the envelope mass fractions are 
qn = 9.5 x 10~ 5 and qn e = 3.1 x 10~ 2 . The core is uni- 
formly mixed with C and O mass fractions of 0.1 and 0.9 
respec tively. Lastly, we use the models of IRenedo et al.l 
(2010), specifically those made publicly available with 
Z=0.001. The envelope thickness of these models varies 
with time and is shown in Figure 4 of IRenedo et al.l 
(2010). The core composition is discussed in section 3.6 
of that same paper with the C/O mass fractions being 
roughly 0.3 and 0.7. 

6.2. Model Comparison 

As can be seen in Figure [T0l the general shape of the 
empirical cooling curve agrees reasonably well with all 
four models. The models all show a linear region with 
a slope of approximately —0.20 followed by a transition 
region, and then a section that agrees with Mestel's value 
of —0.4. From our simple broken power law fit we find 
that the centre of our power law break occurs at an age 
of 104 ± 26 Myrs, and a temperature of 19900 ± 1600 K. 
Fitting the same broken power law to the models yields 
transition ages of 190 Myrs (Fontaine et al.), 129 Myrs 
(Lawlor & MacDonald), 107 Myrs (Hansen), and 250 
Myrs (Renedo et al.). This makes the time of our power 
law transition consistent with both Lawlor & MacDonald 
and Hansen within 2a, but inconsistent with Fontaine 
et al. and Renedo et al. at greater than 3a. How- 
ever, the temperatures at which this transition occurs 
in the models seems to differ from our data. These are 
14,100K (Fontaine et al.), 16,600K (Lawlor & MacDon- 
ald), 16,200K (Hansen), and 13,500K (Renedo et al.). 
These differ at around 2.5c for Lawlor & MacDonald and 



Hansen, and close to 4<r for Fontaine et al. and Renedo 
et al. It is difficult to draw conclusions about how well 
these models fit the data by condensing our entire cooling 
curve down into a few statistics. A more robust method 
involves comparing our observed cumulative distribution 
of temperatures to those that would be predicted by the 
models. In such an analysis, the scale of the x-axis is 
entirely irrelevant since only the relative time spent at 
a given temperature compared to any other temperature 
will determine the shape of the cumulative distribution. 
This approach also compares the distribution over the en- 
tire temperature range, rather than condensing our curve 
down to a few statistics, such as the location of the power 
law break. 

From these cumulative distributions we can calculate 
the Kolmogorov-Smirnov (KS) separation statistic "D" , 
indicating the largest separation between the cumulative 
distribution of our data and that of the model. Next, 
we draw samples of 813 objects from each model distri- 
bution 500,000 times. Here we use a sample size of 813 
rather than our full sample of 887 because only 813 of 
our objects have temperatures that fall between 7,000 K 
and 30,000 K. Each time a random sample of tempera- 
tures is drawn from the model, each temperature is given 
a random shift corresponding to our measured temper- 
ature fit error distribution. All of the temperatures in 
each sample are also given the same shift correspond- 
ing to our systematic temperature error resulting from 
an uncertain distance to the cluster. The uncertainty in 
distance leads to a systematic uncertainty in all of our 
absolute magnitudes and therefore a systematic uncer- 
tainty in all of our fit temperatures. The errors added to 
the model samples correspond to a la distance modulus 
uncertainty of 0.15. This is a conservative estimate, and 
is actually 50% gre ater than the standar d deviation of 
the values quoted in IWoodlev et al.l ()2012[ ) . In each iter- 
ation we calculate the value of D between the sample and 
the model it was drawn from. We then build a distribu- 
tion of these D values and ask how likely it would be to 
draw a sample from these models that produce a value 
of D greater than or equal to that of our data. These 
distributions are shown in Figures [T3] and 1141 

Even with this added uncertainty, our data are ruled 
out as being consistent with any of these cooling models 
at around 3a. Hansen and Renedo et al. are consistent 
at roughly 3a, Lawlor & MacDonald at around 3.5tr, and 
Fontaine et al. at a bit less than 3a. This is due largely to 
the sharpness of the transition in our cooling curve, and 
the temperature at which the transition occurs. None of 
the models show such an abrupt transition to the Mestel 
slope, and so the predicted cumulative distributions dif- 
fer significantly from the observations. It is important 
to note that the age scale determined from the giant 
branch in Section 15.21 can only add further constraint. 
This value will not affect the cumulative distributions 
presented above in any way, and can not make the data 
more consistent with the models. 

We also used othe r cooling models available from 
Ber E'cron etail J20ii|). We carry out the same compar- 
ison described above with the mixed C/O thick atmo- 
sphere models, as well as the mixed C/O thin atmosphere 
models, in addition to the pure C thick atmosphere mod- 
els used above. We found that, given our sample size 
and temperature range, we can not distinguish between 
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Fig. 11. — This plot shows measured cumulative distribution of 
our objects found to have temperatures between 7,000 K and 30,000 
K (813 objects WFC3, 188 objects ACS), as well as those predicted 
by the four models shown in Figure ITDl 
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Fig. 12. — Cumulative distribution of our fit temperatures com- 
pared to the models when assuming values at the extreme ends of 
literature distance measurements for 47 Tuc. 



these models and our data are not consistent with any of 
them. Unfortunately, from this comparison we gain no 
insight into the reasons for the discrepancy between the 
data and the models. It appears this difference cannot 
be easily explained by varying the core composition or H 
atmosphere thickness. 

In addition, we consider the effects of varying the DB 
fraction on the expected cumulative distribution of our 
temperatures and magnitudes (our original approach as- 



sumes our sample is 100% DA stars). We find that we 
can not make the data significantly more consistent with 
the models of iFontaine et al.l (|2001l ) by allowing the DB 
fraction to vary freely. 

A comparison of the temperature distribution from the 
ACS field to each model is shown in Figure Q31 Even 
though the difference between the model distributions 
and ACS data is almost identical to the difference be- 
tween the WFC3 data and the models, due to the smaller 
sample size (188 vs. 813) these data are consistent with 
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13. — Results from our iterative error analysis using data 
from WFC3. The blue line indicates the value of D calculated when 
comparing our data to each model. Our data are inconsistent with 
being drawn from any of these models at the 3<r level. This is true 
even after considering a conservatively large systematic uncertainty 
in distance and a random uncertainty in temperature. 



all of the models at around 2a. 

We also demonstrate in Figure rjj] that the inconsis- 
tency between the models and our data can not be ex- 
plained in terms of a significantly different distance. We 
repeat our fitting procedure while assuming the two most 
extreme dista nce measurements from those summarized 
in Table 1 of iWoodlev etal! (pOlU ). We then build a 
cumulative distribution of fit temperatures and compare 
this to the models. Even centering our distance prior 2a 
from our mean distance of 13.30, our results are incon- 
sistent with any of the models at just below 3a. 

Finally, we perform a simple sanity check to show that 
the difference between the data and the models is not 
a result of bad data in one or more of our filters. We 
calculate the expected cumulative distribution of object 
magnitudes in each of our four filters for each model. 
We then compare this to the observed cumulative dis- 
tributions. These distributions are shown in Figure 1151 
All of these distributions show the data weighted toward 
brighter magnitudes than the models. This agrees with 
the results we find when making the comparison in tem- 
perature space. We also show the comparison to the 
luminosity function in F39QW in Figure [1^1 We also ob- 
serve a discrepancy between the models and data in this 
space as the luminosity function is just the derivative of 
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Fig. 14. — Results from our iterative error analysis using data 
from the ACS field. 

the cumulative distribution (or, more naturally perhaps, 
the cumulative distribution is the integral of the stan- 
dard luminosity function as shown in Equation [5]) In 
this equation the cumulative distribution is calculated 
over the magnitude range [M\, M2] and normalized to 
1. The term $(M) is the luminosity function (dN/dM). 
Calculating the cumulative distribution directly avoids 
the need to bin the data in either magnitude or temper- 
ature. 



jg HM')dM> 



C(M) 



(6) 



7. DISCUSSION 



We have presented an empirical measurement of the 
cooling rate of white dwarfs in the cluster 47 Tucanae. 
Our results suggest two linear regions (in log-space) with 
a transition occuring around 20,000K. In the upper re- 
gion, between approximately 35,000K and 20,000K, the 
temperature of the white dwarfs changes like i~ 019 . This 
slope agrees with all of the models within uncertainties. 
In the cooler region, between approximately 19,000K and 
7,000K, the temperature decreases like t~ 0A2 . This is 
close to Mestel's original value of —0.40 and in fact agrees 
precisely with more thor ough modern c alculations, such 
as those summarized in lHan scn (2004). This region of 
the cooling curve agrees with all of the models as well. 
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Fig. 15. — Cumulative distributions of our object magnitudes in each filter compared with those predicted by the models. 
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Fig. 16. — The luminosity function of our WFC3 data compared 
to the models over the F390W range 21 , 26.5 1. 



However, there is a clear discrepancy between the abrupt 
transition seen in the data and the more gradual transi- 
tion seen in the models. 

Our measurements of both the white dwarf tempera- 
tures and their age on the sequence use a robust maxi- 
mum likelihood fitting method and a thorough treatment 
of all the uncertainties involved. The disagreement of our 
data with the available models, as well as inconsistency 
between the models themselves, suggests differences in 
the input physics involved in the transition region be- 
tween the two power law regimes. The results of our anal- 
ysis can therefore be used to constrain the input physics 



of these models. 

Such cooling models are used to predict the lumi- 
nosity function of white dwarf populations of various 
ages, and are used t o date stellar pop ulations, such 
as the Galactic disk |Oswalt et aLl fl996). old globular 
clusters dHansen et al.ll2007t ). and younger open clusters 
(|von Hippell I2005D . These results therefore have poten- 
tially widescale implications for the age of the Galaxy 
and it's various components. 

Our fit temperatures from all 887 objects from the 
WFC3 data and 292 objects from our ACS data, as well 
as the age for these objects, will be made available online. 

Support for the program GO-11677 was provided by 
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Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under NASA 
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by the NSERC Canada. 
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